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High velocity outflows from young star-forming galaxies in the 
UKIDSS Ultra-Deep Survey 
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ABSTRACT 

We investigate galactic-scale outflows in the redshift range 0.71 $; z ^ 1.63, using 413 
K-band selected galaxies observed in the spectroscopic follow-up of the UKIDSS Ultra- 
Deep Survey (UDSz). The galaxies have an average stellar mass of '-^10^'^ Mq and span 
a wide range in rest-frame colours, representing typical star-forming galaxies at this epoch. 
We stack the spectra by various galaxy properties, including stellar mass, [O ll] equivalent 
width, star-formation rate, specific star-formation rate and rest-frame spectral indices. We 
find that outflows are present in virtually all spectral stacks, with velocities ranging from 100- 
1000 kms^^, indicating that large-scale outflowing winds are a common property at these 
redshifts. The highest velocity outflows (>500 km s~^) are found in galaxies with the highest 
stellar masses and the youngest stellar populations. Our findings suggest that high velocity 
galactic outflows are mostly driven by star-forming processes rather than AGN, with implied 
mass outflow rates comparable to the rates of star formation. Such behaviour is consistent 
with models required to reproduce the high-redshift mass-metallicity relation. 

Key words: galaxies: active - quasars: general - galaxies: evolution 



1 INTRODUCTION 



One of the key unanswered questions in extragalactic astronomy is 
the origin of the galaxy red sequence, and the strong bimodality in 
the co lour-absolute-magnitude plane (e.g. Visvanathan & Sandagd 
1 19770 . This bimodality appears to be in pla ce from the local Uni- 
verse up to red shifts of at least z ~ 2 (e.g . Icirasuolo et al.ll2007l : 



lGeorgakakis"et al. 2008; Chuter et al.ll201l[) . More than half of all 
present-day early-type galaxie s appear to have evo lved into quies- 
cent galaxies by 2 ~ 1 (e.g. lYamada et al J 120051) . From z ~ 1 
to the present day , the red sequence almost doubles in mass (e.g. 
iBundv et alJuOOq) . which is thought to be due to galaxies quench- 
ing their star formation and therefore moving across the colour- 
absolute magnitude plane from the blue cloud to the red sequence, 
from which point on they evolve passively. 

One mechanism for terminating star formation is via galactic- 
scale winds, which eject the gas from a galaxy and deplete the ma- 
terial needed to form stars. Once star formation has been termi- 
nated, the galaxy becomes a passive (quiescent) galaxy, at which 
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point it migrates onto the red sequence. Galactic outflows can 
be driven by stellar processes (e.g. from stella r winds and super- 
novae) or from ac tive galactic nuclei (AGN; e.g. lSilk & ReesI 19981 ; 
iKing et alj|201ll ) and may also play a key role in regulating the 
metal content of a galaxy and enriching the intergalactic medium 
(IGM). 

In starburst galaxies, the energy provided by winds from 
young stars and supemovae is known to drive super-heated bubbles 
of metal enriched plasma (e.g. Heckman et al. 1990). Such bubbles 
expand, sweeping up the interstellar medium as it accelerates to 
velocities of a few hundred km s~^. There is recent evidence (e.g. 
IPiamond-Stanic et alj|2012h to suggest that such winds can reach 
velocities of >1000kms~^ if they have been driven by starburst 
galaxies with a high surface density of star formation. 

In addition to starburst driven winds, feedback from AGN 
is now required in models of galaxy forma tion, to s uppress the 
formation of the most massive galaxies (e.g. iBenson et al. 200^. 
A number of models use AGN feedback in the form of galactic- 
scale outflows, sufficient to remove the majority of gas, quench 
star formation and thereby create massive , passive galaxies (e.g. 
iDi Matteo et alj|2005l ; [Hopkins et alj|2005h . Low-luminosity AGN 
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activity may then suppress any further star formation as it ' sweeps 
up' the remaining gas which would o therwise have been used 
to form more stars (e.g. iHopkinsI 120111) . Alternative models use 
'radio-mode' feedback in the form of expa nding bubbles within the 
hot intergalactic medium (e.g. lBower et al . 2006). The precise form 
of AGN feedback remains unclear, and a combination of radio- 
mode feedback and AGN- driven winds may be required to suppress 
star formation sufficiently jvan de Voort et al.ll201 II) . 

Galactic winds can be detected through the study of many dif- 
ferent emission and absorption lines, but in particular, spectra in the 
restframe UV are ideal for probing the interstellar medium (ISM). 
This part of the spectrum is rich in interstellar absorption lines, 
such as Mgll (2796, 2803 A) and Fell (2344, 2374, 2382, 2586, 
2600 A), which can be used to trace cool interstellar gas. Outflows 
with velocities of up to ~600kms^^ are often det ected in star- 
burst galaxies and appear to be present at all epochs (IWeiner et al.l 
[2009; Rubin et al. 2010; Coil et al. 2011) and are thought to be the 
product of supernova explosions. Outflow signatures can also be 
detected i n the far UV with absorption lines such as C IV, Sill, 
Si IV etc ISteidel et al.ll2004). Outflows with higher velocities are 
thought to be AG N-driven (Thacke r etai1l200a ISato et al.ll2009l ; 
ISturm et al.ll201 II) . although (as previously menti oned) there is re- 
cent evidence from lOiamond-Stanic et al.l ( 120120 to suggest they 
could also be starburst driven. 

iTremonti et al.l ( 120071) discovered relic outflows with veloci- 
ties > 1 000 km s ~ ^ in 1 out of 1 4 high resolution spectra of highly- 
luminous post-starburst galaxies at a redshift of z ~ 0.6 using the 
Mg II AA 2795.53, 2802.71 A doublet. Post-starburst (also referred 
to as E+A) galaxies are those in which star formation has recently 
been terminated and are assumed to be in the transition to becom- 
ing an early-type galaxy. The high occurrence of high-velocity out- 
flows in bright, post-starburst galaxies suggests that these outflows 
could be t he quenching mechanis m responsible for terminating star 
formation. ITremonti et al.l ( 12007 ^ suggest that these winds may be 
AGN-driven, given the high velocity of the outflows. It is unclear, 
however, if this is the mechanism by which more typical galaxies 
migrate to the red sequence. 

Work on winds and galactic-scale outflows has a lso been car- 
ried o ut at higher redshifts (z > 2). For example, IPettini et al.l 
( l2002h studied galactic winds by observing a gravitationally lensed 
galaxy with a redshift z = 2.72. They utilised UV absorption lines 
to measure outflows and found that although the bulk of the inter- 
stellar medium has a velocity ~ 250 km s~ ^ , there are also outflows 
present with velocities of at least three times this value. Pettini et 
al. deduced from this work that (if the outflow is spherical in geom- 
etry) the mass outflow rate is greater than the star-formation rate. 
However, they could not determine whether the gas stays gravita- 
tionally bound to the g alaxy and henc e whet her the winds could ter- 
minate star formation. IShaplev et al.l (l2003h investigated the pres- 
ence of outflows in Lyman break galaxies (LBGs) at z '--^ 3 by com- 
paring redshifts derived from interstellar lines with Lyman-alpha 
emission. They found outflows exist at this redshift with a wide 
range of velocities and the average offset between the velocities of 
the interstellar absorption lines an d the Lyman-alpha emi ssion was 
found to be ~650kms~^. Finallv. lSwinbank et al.l j2005l) observed 
a SCUBA galaxy at 2 = 2.4 and found a gas outflow with a veloc- 
ity of 220kms~^, which they interpret as a starburst-driven wind. 
Outflows at higher redshifts seem extremely common, but outflows 
in galaxies with an intermediate redshift have been studied veiy 
little until recently. 

In this work we investigate galactic scale outflows up to red- 
shifts z ~ 1.6, the epoch when the red sequence is still being 



rapidly assembled and AGN activity peaks. We use a sample of 
A'-band selected galaxi es, probing more typical galaxies than the 
work byiTremonti et al. (2007). Our work differs from other recent 
studies ( Weiner et al. 2009; Rubin et al. 2010; Coil et al. 201 1) by 
the inclusion of many very red galaxies at z > 1, including many 
on the galaxy red sequence. We use 413 galaxy spectra with red- 
shifts in the range 0.71 < z < 1.63, stacking the spectra by galaxy 
colours, stellar mass and spectral properties. Our aim is to deter- 
mine the relationship between galaxies properties and the occur- 
rence of high-velocity outflows, as measured from the Mg II ab- 
sorption feature. 

Where relevant, we adopt a concordance cosmology in our 
analysis; Q.m = 0.3, Qa = 0.7, h = Ho/100 kms'^Mpc"^ = 
0.7 and erg = 0.9. Section 2 describes the photometric and spectro- 
scopic data used in this work, while Section 3 describes the method- 
ology of our stacking and line-fitting analysis. Section 4 presents 
the results of our work based on an analysis of Mg II feature, while 
Section 5 discusses the implications for terminating star formation. 
Section 6 presents a summary of our conclusions. 

We separate the spectroscopic data used in this analysis into 
those obtained with the FORS2 instrument and those obtained with 
VIMOS, due to the differing resolutions of the instruments. We de- 
scribe the VIMOS methodology in the main body of the text and 
discuss the details of the FORS2 analysis in the appendix. 



2 SAMPLE SELECTION 

2.1 UDS Photometric Data 

To investigate the prevalence and velocity structure of outflows, our 
galaxy sample is drawn from the Ultra Deep Survey (UDS). The 
UDS is the deepest of the five surveys that comprise the UKIRT 
Infrared Deep Sky Survey (UKIDSS; Lawrence et al. 2007). The 
UDS covers an area of 0.77 square degrees, centred on the Subaru- 
XMM Deep Field (SXDF), and is the deepest near-infrared survey 
over such a large area to date. In this work we use data from the 
UKIDSS data release 8 (DR8), reaching depths of K = 24.6, H = 
24.2, J = 24.9 (AB, 5a, 2" aperture). In addition to near infrared 
data, the field is covered by deep optical data in the B, V, R, i' and 
z'-bands from the Subaru-XMM Deep Survey, achieving depths of 
B = 28.4, V = 27 .8, R = 27.7, i' = 27.7 and z' = 26.7 
dFurusawa et al.l2008l . AB, 3a, 2"). Spitzer data reaching 5a depths 
of 24.2 and 24.0 (AB) at 3.6/^m and 4.5/im respectively is pro- 
vided by the Spitzer UDS Legacy Program (SpUDS, PLDunlop) 
with U-band data taken with the Canada-France-Hawaii Telescope 
(CFHT) using the Megacam instrument (Uab = 25.5; Foucaud et 
al. in prep). Although spectroscopy is used for the majority of the 
analysis, we also utilise UDS imaging to calculate galaxy stellar 
masses and to compare star-formation rates derived from emission 
lines with SED fitting. 



2.2 UDS Spectroscopic Data 

The spectra used in this analysis were taken from the ESO Large 
Programme 180.A-0776 (UDSz; PI: Almaini), using a combina- 
tion of the VIMOS and F0RS2 instruments (Visible Multi-Object 
Spectrograph and FOcal Reducer and low dispersion Spectrograph 
respectively), installed on the ESO VLT. 

The UDSz programme was designed to target a large K- 
band selected sample of ~3500 galaxies across the UDS field, 
with photometric redshifts used to preferentially target galaxies at 
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Zphot > 1- VIMOS was used to target optically brighter systems 
(i'AB < 24 or Vab < 25) whilst F0RS2 was used for fainter, red- 
der systems (iab < 25 and Vab > 24) with deliberate overlap in 
the selection criteria. 

The UDSz survey used eight pointings of VIMOS in multi- 
object mode to sample an area of 0.5 sq degrees across the UDS 
field, typically observing over 300 targets at each position. For each 
of these pointings we obtained 4.5 hours exposure with the LR Blue 
grism and 2.6 hours with the LR Red grism with resolving powers 
of R — X/5X — 180, 210 respectively, with the majority of tar- 
gets observed with both grisms to maximise wavelength coverage. 
With the LR Red and LR Blue grisms we obtained an average RMS 
of '^ 1.1 A and ~0.93 A in the wavelength calibration respectively, 
and for the FORS2 data, this decreases to ~0.24 A . There were 
20 FORS2 fields coincident with the central part of each VIMOS 
tile; each of which received 5.5 hours exposure with the 3001 grism, 
typically observing 35-40 faint galaxies per field with a resolving 
power of R = 660. In total, 2881 objects were targeted with VI- 
MOS and 802 with FORS2. 

The VIMOS and FORS2 spectra were reduced separately 
at different i nstitutions. For more inf ormation on the F0RS2 re- 
duction, see ICurtis-Lake et al.l ( l2012h . The VIMOS spectra were 
firstly extracted and reduced using the VIMOS Inte ractive Pipeline 
Graphical Interface (VIPGL lScodeggio et al.l2005h . The zeroth or- 
der lines and known sky lines were then removed from the ID spec- 
tra. The spectra were trimmed to remove regions badly affected 
by fringing and then VIMOS LR Red and LR Blue spectra were 
spliced together (where applicable) to yield typical spectra with a 
total wavelength range of 3700-9500 A. The spectra were resam- 
pled and spliced together using the mean S/N in the overlap region 
(6000-6600 A) to weight the spectra. The average continuum S/N 
measurement was found to be 6.46 A ~^ for the red spectra and 4.49 
A ^^ in the blue. The spectra we use in this analysis are those with 
highly secure redshifts (see below), for which the S/N is typically 
higher. In order to obtain the high S/N required for the detailed 
study of the Mg II absorption line, however, we create composite 
spectra by stacking the individual galaxy spectra in the restframe 
(see Section 3.2). 

2.2.1 Spectroscopic Redshifts 

Redshifts for t he spectra were obt ained using two packages: SG- 
NAPS and EZ JFumana et al.ll2008r) . EZ determines redshifts using 
the cross-correlation of templates, the solutions of which can then 
be checked manually. We used a selection of default templates and 
also our own templates which were constructed after a first pass of 
redshift determination. These templates will be presented in Brad- 
shaw et al. (in preparation). Overall, we obtained 1512 highly se- 
cure redshifts from the primary iC'-selected sample spanning a red- 
shift range < z < 4.8, including several hundred in the 'redshift 
desert' (1.4 < z < 2.2). A number of z > 6 candidates were also 
confirmed (see Curtis-Lake et al. 2012). 

For the outflow analysis it was important to use only those 
spectra with a secure, accurately-determined redshift, and there- 
fore we chose to use only spectra with strong [O II] emission at 
3121.5k. , which is easily identifiable and close to Mg II in wave- 
length. We also further limit the sample by selecting only those 
spectra with [O II] of sufficient strength to yield a redshift accu- 
rate to ~I00kms^^, as determined by fitting to the centroid of 
the emission line. An eiTor of this magnitude is typically twice the 
error in the VIMOS wavelength calibration. By using [O II] as a 
zero-point to define the systemic redshift we ensure that the posi- 



tion of Mg II is measured consistently for all spectra. Although we 
cannot resolve the lines, [O II] is actually a doublet with emission 
lines at 3726. 1 A and 3728.8 A . NV e assume a doublet ratio of 8:9, 
as observed bv lWeiner et al.l (120091) using higher-resolution spectra 
for a sample of galaxies with similar redshifts and stellar masses. 
This is equivalent to an electron density of ~5 x 10* m~"^. The cen- 
tre of the [O II] emission therefore lies at 3727.5 A . We use this 
wavelength throughout the rest of our analysis. 

In order for both [O II] and Mg II to be visible in the ob- 
served spectra we effectively restrict our analysis to redshifts in 
the range 0.71 < z < 1.41 for the VIMOS data, extending 
up to z ~ 1.63 for the FORS2 data. Finally, the spectra identi- 
fied as clear AGN were removed from thi s data analysis; thi s in- 



cluded those i dentified as X-ray AG N by lUeda et al.l (|2008|) and 
radio AGN by ISimpson et al] ( |2006|) , (see iBradshaw et al.ll20lll 
for more information) as well as those spectra displaying emis- 
sion lines commonly identified as originating from AGN, such as 
[C IV] and C III] emission. Spectra with badly subtracted skylines 
were also excluded from this analysis, leaving a sample total of 413 
galaxy spectra, with 341 from VIMOS and 72 from FORS2. The 
413 galaxies in the sample have an average stellar mass of ~10^'^ 
M0 and display a wide range of rest-frame colours (see Section 
4). Given the primary selection in the A'-band, our sample will be 
broadly representative of typical galaxies at this epoch. The most 
passive galaxies (with little or no star formation) will be somewhat 
under-represented, due to the requirement for [O II] emission in our 
sample selection, but we note that the F0RS2 sample in particular 
contains many galaxies located on the galaxy red sequence (see 
SectionHjl. 



2.2.2 Stellar Masses 

The stellar masses and magnitudes used in this work are measured 
using a multi-colour stellar population fitting technique where we 
fit to the UBVRizJHK bands and IRAC Channel 1 and 2 bands. 
A large grid of synthetic SEDs are constructed from the stellar pop- 
ulation models of Bruzual & Chariot (2003), assuming a Chabrier 
initial mass function. The star-formation history is characterised by 
an exponentially declining model with various ages, metallicities 
and dust extinctions. These models are parametrised by an age of 
the onset of star formation, and by an e-folding time such that 



SFR{t)r^ SFRo xe--. 



(1) 



where the values of r can range between 0.01 and 13.7 Gyr, while 
the age of the onset of star formation ranges from 0.001 to 13.7Gyr. 
The metallicity ranges from 0.0001 to 0.1, and the dust content is 
parametrised by r^, the effective 1/-band optical depth for which 
we use values -r„ = 0, 0.2, 0.4, 0.6, 0.8, 1.0, 1.33, 1.66, 2, 2.5. 

To fit the SEDs we first scale them to the A'-band magnitude 
of the galaxy we are fitting to. We then fit each scaled model tem- 
plate in the grid of SEDs to the measured photometry of the galaxy. 
We compute the chi-squared values for each template and select the 
best fit. From this we obtain a best-fit stellar mass and best-fit ab- 
solute magnitud es. Further details of the stellar mass determination 
can be found in lHartlev et al.l ( 120130 and Mortlock et al. (submit- 
ted). 



3 CREATING COMPOSITE SPECTRA 

We use the Mg II absorption feature as a tracer of outflowing galac- 
tic winds because it is a sensitive tracer of cool gas in the ISM. 
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Figure 1. The co-added stack of the 341 VIMOS spectra used in this work, all selected to ensure that [O II] and Mg II were within the observed spectral range. 
Other emission and absorption features are also shown. Flux is shown in ergs s~^ cm~^. 



To study the properties of the Mg II AA 2796, 2803 doublet, we 
performed a stacking analysis with our spectra, as there are only a 
handful of galaxies which are bright enough to individually mea- 
sure the structure of the Mg II profile. Higher S/N composites allow 
us to study the relationship between the average outflow properties 
and a range of galaxy characteristics. 

The stacking analysis was conducted for the VIMOS and 
F0RS2 data separately as the resolutions of the two instruments 
are different. For clarity, we place details of the FORS2 analysis in 
the appendix and describe here the main analysis with the VIMOS 
spectra, as these data comprise the bulk of the sample (341 VIMOS 
spectra compared to 72 F0RS2 spectra). 



3.1 Co-addition of UDSz Spectra 

The redshifts obtained during the data reduction process were used 
to shift the spectra to the restframe whilst oversampling the spectra 
onto a finer dispersion axis so as not to lose information. In order 
to confirm the reliability of the zero-point redshift derived from the 
[On] line at 3727.5 A, we calculate the offset between [Oil] and 
the Call K absorption line at 3933.7 A , which should be unaffected 
by any systemic redshift offset from outflowing gas, as it is primar- 
ily a stellar absorption line. From the composite spectrum of the 
VIMOS sample we find that the centroid of the Call absorption line 
is offset by 0.2 A (~ 20 km s^^) from the [O II] emission, which is 
well below the value of our wavelength calibration error. We also 
compare the velocity offset between [O II] and C III], as this feature 
is towards the bluest end of our spectra and hence provides a good 
test of our calibration uncertainties over a wide wavelength range. 
The Cm] A 1909 A emission is known to be a doublet (AA 1906.1, 
1908.7 A), where the flux density of the line s depend particularl y 
on the electron density of the emitting gas JKeenan et alj|l992h . 
With our data, we cannot resolve the C III] lines, so we use a ratio 



of 2:3 with a centre at 1907.7 A , corresponding to an electron den- 
sity of ~5x 10* m""* within star-forming regions, consistent with 
the elec tron density we implicitly assume for the [Oil] doublet ra- 
tio (e.g. lTalia et al.l2012r) . The C III] emission, was found to have an 
offset of 0.3 A (~ 30km s~^) from [O II], which is again below the 
error in our wavelength calibration. Similar values were obtained 
using the other composite spectra used in this work. These checks 
suggest that the [O II] line provides a reasonable measure of the 
systemic redshifts. 

Two subtly different routines were used to construct the com- 
posite spectra. In both methods, the continuum levels of all spectra 
were scaled by fitting the continuum in the restframe wavelength 
range 2850-3050 A , a featureless window near Mg II (2800 A ). 
The spectra were co-added, firstly by taking the mean pixel value 
at each wavelength increment (thus ensuring a high S/N) and sec- 
ondly by using the median pixel value at each wavelength incre- 
ment (thus ensuring the signal from any unusual galaxies do not 
dominate the co-added spectrum). The composite spectrum of all 
341 VIMOS galaxies in our sample (using the median method) is 
shown in Figure [T] We find that the chosen method of stacking has 
no significant impact on the results presented, so we use the median 
stacking method throughout this work. 

A close-up of the Mg II absorption profile from the stack of 
all 341 VIMOS galaxy spectra is shown in Figure |2] By fitting a 
Gaussian curve to the profile with its centroid at km s^^ we find 
clear evidence for excess absorption blueward of the fitted line, sug- 
gesting a component of outflowing gas. By fitting a Gaussian to 
the whole Mgll profile (without fixing the centroid at 0km s~^), 
we measure a blueshift of ~I.5A which corresponds to an out- 
flow velocity of ~150kms^^. However, this simple fit is likely to 
underestimate the velocity of any outflowing component if there 
is a significant component of non-outflowing gas in the stack. We 
therefore use a method which fits two Gaussian profiles (as detailed 
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Figure 2. The Mg II absorption profile from a co-added spectrum containing 
all 341 VIMOS galaxies in our sample. The red curve shows the best-fitting 
Gaussian curve, centred at kms~^. Clear absorption blueward of the line 
is visible. By fitting a two-component model (see Section 3.2.1), we obtain 
a velocity of ~340kms~^ for the outflowing component. 



in Section 3.2) to parameterise the intrinsic and outflowing com- 
ponents of Mg II. This method yields an average wind velocity of 
340 km s^ ^ for the outflowing component. 



3.2 Isolating Outflowing Mgll 

In this Section we describe the decomposition of the Mg 11 absorp- 
tion profile into two separate components; a non-outflowing com- 
ponent (at the systemic velocity) and an outflowing component, 
with the velocity determined from a fitting procedure. 

We assume that, in a composite spectrum of galaxies, there 
will be a sufficiently large number of spectra such that the absorp- 
tion due to non-outflowing galactic gas will have a mean velocity 
of zero, matching the systemic velocity determined from the [O II] 
line. On the other hand, the outflowing gas from a galaxy will al- 
ways have a zero or negative velocity, depending on the geometry 
of the outflow. We would not be able to see the positive velocity 
component of the outflow (should the outflow be symmetrical) due 
to the lack of continuum light sources along the line of sight. 

It should be noted that the Mg II equivalent widths that we de- 
rive below from the composite spectra do not necessarily provide a 
means for determining the column density of the absorbing gas. Al- 
though it is not obvious from our low-resolution spectra, the Mg II 
absorption in an individual galaxy is typically very close to satu- 
rated (Weiner et al. 2009). Therefore, in this work, the absorption 
depths and measured equivalent width provide only an approximate 
indicator of the fraction of galaxies at a given velocity in the com- 
posite spectrum. 



3.2.1 Methodology: Extracting the Intrinsic and Outflowing 
Mgll Components 

Firstly, to prepare the spectra for fitting, they were trimmed to a 
wavelength region that includes Mgll and [Oil] (2550-3800 A). 
The continuum in this region was fitted with a spline function and 
removed by dividing the spectrum by the fit. 

To isolate and examine the intrinsic and outflowing gas, we 
complete an analysis which involves a careful fitting of two model 



Mgll Gaussian absorption components; one representing the in- 
trinsic (non-outflowing) component of Mg II (at zero velocity) and 
the other representing the outflowing component (at negative ve- 
locities). For the intrinsic Mg II, we centre a Gaussian curve at 
Okms~^, all owing only the line depth to vary. A similar method 
was used by I Weiner et al.l ( 120091) to successfully isolate the sys- 
temic component from higher resolution spectra. The depth and 
velocity of the outflowing component were free parameters. 

Both Gaussians were restricted to have the same FWHM, to 
match the measured instrumental resolution of VIMOS. This width 
was determined using the [O II] emission line in the composite 
spectra at 3727 A, allowing for the doublet nature of this feature. 
As an alternative method we also determined the FWHM by allow- 
ing a free fit to the red side of the absoiption profile. Both methods 
gave consistent results, as noted below. 

The two components are fitted jointly (using a least-squares 
fit) to the absorption profile of Mg II. The results of this analysis 
can be seen (for example) in the right-hand column of Figure |4] 
where the red curve represents the intrinsic Mg II, the green curve 
represents the outflowing Mg II, and the blue curve shows the over- 
all fit to the absorption profile. 

In addition, to test these results, we also created a more com- 
plex model of the absorption profile, modelling the Mg II doublet 
with two components centred at 2796 A and 2803 A, each with 
a = 60kms~^ and using a doublet ratio of 1.1:1, as observed in 
higher resolution spectra of galaxies at z '^ 1 by Weiner et al. 
(2009). This in turn was convolved with another Gaussian to model 
the response function of the instrument. The resulting model was 
then fitted to the absorption profiles with no noticeable difference 
in results between this and the two-Gaussian method previously 
described. 

In order to test the robustness of this method further, we also 
determined the instrumental FWHM by first fitting a Gaussian pro- 
file to the red side of Mg II, but otherwise following the same pro- 
cedure as detailed above. This was in order to test that the FWHM 
of [O II] was a fair representation of the instrumental resolution (af- 
ter accounting for its doublet nature). This alternative Mg II profile 
was then fitted to the data in Section 4. No significant differences 
were found between the two methods. The maximum variation of 
the outflowing velocity component from the method using a fixed 
FWHM is llOkms^^ and the average deviation between the two 
methods is 60 km s~^, where the Mg II with the free-fitting FWHM 
extracts higher velocity outflows (in the vast majority of cases). In 
every case, the two methods are both consistent with each other 
within the margin of error calculated using the method in Sec- 
tion [33] For the remainder of this analysis we therefore present 
results determined by fixing the FWHM at the instrumental value 
determined from the [O II] emission line, as described above. 



3.2.2 Emission Filling and Resonant Scattering 

In the work by Weine r et al.l ( 12009 !), it was noticed that there was an 
emission feature present in the spectra of low-mass galaxies with 
bluer UV slopes on the red side of Mg II. The emission was at- 
tributed to possible AGN activity and all galaxies showing signs 
of this feature were removed from subsequent analysis. More re- 
cent papers on the topic have discovered that the emission may 
be caused by resonance-line scattering from receding gas on the 
far side of the galaxies that host outflows, resulting in a P-Cygni- 
like line profile with redshifted emission and blueshifted absorp- 
tion (e.g.lProchaska et al. l l20I ll ; iRubin et aklboill ; lErb et al.ll2012l ; 
iKomei et alj|2012l) . This phenomenon can lead to significant emis- 
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flow velocities (unless the emission component is extreme). At ve- 
locities below 300kms~^, however, the outflows we infer may, in 
principle, be artificially boosted to higher values. 

Finally, we compared the results of our simulated spectra with 
the real spectral stacks to determine if a strong emission component 
would be visible even at the low resolution of our VIMOS spectra, 
allowing for the noise in the data. We concluded that a strong Mg II 
emission component (> 35% of the absorption strength) would 
always be clearly visible given our signal to noise. The clearest ev- 
idence for a strong emission component was found in the low-mass 
stack presented in the top panel of Figure [7] Comparison with our 
model spectra suggest that the emission feature has an equivalent 
width comparable to ~65% of the absorption component. Based 
on Figure |3] therefore, we believe the outflow velocity can be ex- 
tracted down to an approximate limit of ~300kms~^, which is 
similar to the outflow velocity we extract for this spectrum. As our 
other composite spectra show little or no signs of emission filling, 
we determine that our method is generally sufficient for the rest of 
our analysis. 



Figure 3. This figure shows the results of simulations to quantify the impact 
of a P-Cygni-like emission feature on the Mg II profile, which may affect 
our ability to determine an accurate outflow velocity. A range of outflow 
velocities and emission strengths are modelled. The percentages refer to the 
strength of the emission compared to the absorption. We find that outflow 
velocities can be re-extracted for low emission strengths, but velocities are 
systematically boosted upwards when the emission feature is very strong. 



sion underlying the absoiption, which can in turn make the outflow 
velocity hard to determine. Here we investigate the effect that emis- 
sion filling has on our ability to extract an accurate outflow velocity 
from our data. 

We obtain a general insight as to the effects of emission filling 
by simulating artificial spectra of Mg II with a P-Cygni-like profile 
whilst varying both the amount of emission and the outflow ve- 
locity. We then attempt to re-extract the outflow velocity using the 
Mg II model used in our analysis (as described in Section 13.2. U . 
We begin by mode lling the P-Cygni profile shown in Figure 1 of 
iRubin et al.l (1 20111 ) with four overlapping Gaussians; two repre- 
senting emission and two representing absorption (due to the Mg II 
doublet). The profile in this paper shows extremely strong emis- 
sion at a ratio of nearly 1 : 1 between the emission trough and the 
absorption peak. We then convolve this model with a Gaussian to 
match the instrumental resolution of our VIM OS spectra, a l lowin g 
for the intrinsic resolution of the spectrum in iRubin et al.l (120111) . 
The strength of the emission components are progressively reduced 
until the spectrum consists of a pure absorption profile, and we also 
vary the velocity of the outflow from lOOkms"^ to 1000 km s~^. 

Figure |3] shows the velocities that we re-extract from the arti- 
ficial spectra for increasing absorption/emission ratios and increas- 
ing outflow velocities. For spectra in which the emission is 20% the 
strength of the absorption (or lower), we can easily extract the same 
outflow velocities that are input to the model. For spectra with an 
emission strength varying between 20% and 70% of the absorption, 
we must utilise caution in our analysis of any velocities lower than 
~300kms~^. Finally, for emission which is the same strength as 
the absorption, we cannot accurately re-extract any velocity below 
800kms"\ 

Inspection of our fits suggest that the outflow velocity we ex- 
tract relies heavily on the blue side of the Mg II profile, and hence 
the results are largely unaffected by emission filling at high out- 



3.3 Errors 

In order to obtain an estimate of the error involved in extracting 
a ineasurement for the equivalent width and velocity of the out- 
flowing gas, we undertake a bootstrap analysis of our data. For ev- 
ery composite spectrum, we make an additional 100 spectra from 
sub-samples which are created by choosing a random selection of 
spectra within the same subset, but with replacement. All spectra 
are then resampled onto the same wavelength axis and the standard 
deviation on each pixel is calculated. The fitting procedure from 
Section |3.2.1| is then repeated for each of the bootstrapped spectra, 
and we take the error of the velocity and column density to be one 
standard deviation from the mean of the 100 bootstrapped samples. 



3.4 Single-Gaussian Fit 

In addition to the double-Gaussian fits described above, throughout 
this paper we also present the results of single-Gaussian fits, fixed 
at the systemic redshift and fitting only to the red side of the line. 
These fits allow the velocity and extent of any outflowing compo- 
nent to be visualised from the residuals blueward of the fitted line. 
We follow the procedure in Section |3. 2. ll to fit a Gaussian profile 
which is centred at kms~^ using the instrumental FWHM de- 
termined from the [O II] emission line. The depth of the profile is 
allowed to vary. The single-Gaussian fits help to illustrate whether 
an additional blue-shifted component is required. 



4 OUTFLOWS IN GALAXIES STACKED BY PROPERTY 

In this section, we discuss the properties used to categorise the 
UDSz galaxies into stacks and the results of our subsequent fit- 
ting analysis. This section contains results from both the VIMOS 
sample and the F0RS2 sample, although the detailed methodology 
for analysing the F0RS2 spectra is contained in the appendix. 



4.1 Star Formation 

The equivalent width of [O II] emission at 3727A has long been 
recognised as an approximate indicator for the strength of s tar for- 
mation (e.g. lRosa-Gonzalez et alj2002l : lKewlev et alj2004l ). In the 
absence of significant dust or AGN activity, strong [O II] emission 
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Figure 4. The left panels show stacked spectra, arranged by decreasing [0 11] equivalent width, focussing on the wavelength range 2600-3800 A . The central- 
left panels show a close-up of the Mg II absorption feature for each stack, along with the best-fitting single Gaussian (in red), obtained by fitting to the red 
side of the line only and fixed at zero velocity. The central-right panel shows the errors associated with each absorption profile, which were calculated using 
a bootstrap method. The right-hand panels show the breakdown of the absorption line fits used to calculate the velocity of the outflowing Mg II, based on a 
two-component model. The red curve is the non-outflowing component, the green curve is the outflowing component, and the blue line represents the best-fit 
to the Mgll line (the sum of the stationary and outflowing components). Flux is shown in ergs s^^ cm^^. 
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Figure 5. The velocity (left panel) and the equivalent width (light panel) obtained for the outflowing component of Mg 11 from the two-component flts to the 
stacked spectra shown in Figure 3, displayed as a function of [O II] EW. The green points are based on the FORS2 data and the red points are based on the 
VIMOS sample. Data are plotted at the mean value of the [O II] EW within a stack. 
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indicates a significant burst of recent star formation. Tiierefore, in 
order to investigate tlie relationsliip between star formation activ- 
ity and Mg II outflows, we measured the equivalent width of [O II] 
emission at 3727A for all VIMOS spectra. The spectra were then 
divided into four sub-samples; those showing the most active star 
formation (EW=30-50 A), two intermediate samples (EW=20- 
30 A) and EW= 10-20 A) and finally those with weak star for- 
mation (EW< 10 A ). These bins were chosen to ensure a sufficient 
number of spectra in each stack (in this case there are 37, 5 1 , 84 and 
88 spectra in each stack respectively) whilst encompassing a large 
range of equivalent widths. The co-added spectra can be seen in 
in the left-hand panels of Figure |4] A single-Gaussian fit was then 
performed on the Mg II profile in each stack, centred at km s~^ 
(representing the non-outflowing component of the Mg II), where 
the curve was fitted only to the red side of the Mg II absorption line 
(see Section [J!4t . The excess Mgll and hence any outflowing ma- 
terial can then clearly be seen from the absorption blueward of this 
profile. 

In order to illustrate the errors involved in the fitting process, 
the central-right-hand panels of Figure |4] show the errors on the 
stacked spectra through a bootstrap analysis, which was detailed in 
Section [33] The errors are noticeably larger in the stacks containing 
fewer spectra. 

The velocity and equivalent widths of the intrinsic and out- 
flowing components of Mg II in each sub-sample are then derived 
from the method in Section l3.2.1l using two-Gaussian fits, and are 
shown in the right-hand panel. The equivalent width and the ve- 
locities of the Gaussians representing the outflowing component 
were extracted and plotted in Figure|5]in red, along with equivalent 
points extracted separately from the F0RS2 data (in green). The 
errors are calculated from the methodology in Section [33] Those 
spectra with a high [O II] equivalent width clearly show the highest 
velocity Mg II outflows, but we find no correlation between [O II] 
equivalent width and the equivalent width of the outflowing com- 
ponent of Mgll. 



4.2 Mb Magnitude 

We divide the spectra into four sub-samples dependant on their ab- 
solute B-band magnitude; a bright sub-sample with Mb < —21, 
two intermediate sub-samples —21 ^ Mb < —20.6 and —20.6 ^ 
Mb < -20.2, and those spectra with A/s ^ -20.2. Each of these 
four samples were then co-added (see methodology in Section UlT ). 
and curves were fitted in order to extract the intrinsic and outflow- 
ing Mgll (see Section 3.2.1). For the sub-sample with the bright- 
est B-band magnitude, an outflow velocity of 475 ± 114kms~^ 
was extracted. Wind velocities of 350 ± 73kms~^ and 389 ± 
23 km s~^ were obtained for the two intermediate samples. We do 
not extract a velocity for the faintest sub-sample, as the spectrum 
is too noisy to fit to. We plot these results on a graph adapted from 
iTremonti et al.l (2007). which compares the velocity shift to the av- 
erage B-band magnitude of the galaxies (Figure |6j. Our findings 
compare well with the findings of previous authors, with outflow 
velocities at the upper end compared to local starburst galaxies and 
ultra luminous infrared galaxies (ULIRGS). 



4.3 Stellar Mass 




In thi s section we investigate the claim made by IWeiner et al.l 
( 120091) for a positive relation between galaxy stellar mass and out- 
flow velocity. In order to match our sample to their work, we first 



Figure 6. This figure, adapted from Tremonti et al. (2007), shows the rela- 
tionship between the velocity offset of low ionisation absorption lines and 
the B-band absolute magnitude. Our stacked data are represented by the 
black squai'es and correspond to the top end of the velocity offsets found in 
local galaxies classified as starbursts. 



limit our sample to blue galaxies with re stframe (U — B) < 0.5 
to broadly match the DEEP2 sample from lWeiner et al.l (120091) . We 
cannot perform a similar analysis on the FORS2 data, as there are 
only 27 spectra with a {U — B) < 0.5, which are insufficient for a 
robust fit. 

We divide our sample into three sub-samples, those with 
a stellar mass below lO^'^^'M©, those with a mass higher than 
IO^'^^Mq, and those with a mass between the two. We then cre- 
ate three composite spectra from these sub-samples and follow the 
same methodology as before to measure the velocity and equivalent 
width of outflowing material. As illustrated in Figure|7] we find ev- 
idence for higher velocity outflows in the g alaxies with a highe r 
mass. This trend is broadly consistent with IWeiner et al.l (|2009|) . 
despite the different methods used to extract the velocities of the 
outflows. 

We then put the red galaxies back into the sample and co-add 
the galaxies in the same three sub-samples as above. We follow the 
same methodology used throughout this paper to derive the veloc- 
ity of the outflowing gas. The results are then also plotted on Fig- 
ureQin red. There are no red galaxies in the lowest mass sample 
and therefore the two points representing the co-added spectrum 
with the lowest mass are the same. We see a significant drop in the 
velocity of the outflow when we include the red galaxies, which 
suggests it could be colour and star-formation which most strongly 
determines the wind velocity. 

As mentioned in Section |3.2.2| we are cautious about extract- 
ing velocities using our two-component method for those spectra 
with signs of Mg II emission fil ling. This is a comm on phenomenon 
for galaxies of lower mass (e.g. lWeiner et al . 120091 : lErb et al.l2012h 
and we do see emission in Figure |7] slightly redward of the Mgll 
absorption in our composite spectrum of low mass galaxies. As pre- 
viously discussed, we model the Mg II emission and absorption to 
find that the emission is ~65% the strength of the absorption com- 
ponent and therefore we must use caution in the interpretation of re- 
sults if the extracted velocities are below ~300km s~^ (Figure|3]l. 
In this case, the extracted velocities from the spectra composed 
from higher mass galaxies are above 700km s^^, which suggests 
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Figure 7. The left-hand panels show a close-up of the Mg II absorption feature for three spectral composites arranged by stellar mass, along with the best-fitting 
single Gaussian (in red), which represents the restframe Mg II. The central panels show the breakdown of the absoiption line fits used to calculate the velocity 
of the outflowing Mg II. The red curve is the non-outflowing component, the green curve is the outflowing component, and the blue curve represents the best-fit 
to the Mgll line (the sum of the stationary and outflowing components). These composite spectra are based on relatively blue galaxies with (U — B) < 0.5 
to allow a comparison with similar analyses in the literature (see Section 4.3). The right panel shows the equivalent width and velocity of the outflowing 
components as a function of stellar mass. In red we plot the complete VIMOS sainple and in blue we plot blue galaxies with (U — B) < 0.5 only. The mean 
stellar mass in the lowest-mass and medium-mass bins are the same, therefore for clarity we offset the points from each other 



that there is a definite trend with mass, even if the extracted veloc- 
ity of the lowest mass composite is potentially inaccurate. 



4.4 Star-formation Rate 

Here we use [O II] luminosity as a proxy for star-formation rate to 
determine whether there is a relationship with outflow velocity. We 
calculate the star-form ation rate using the following equation from 
lGallagheret"ZI ( ll989l) : 

SFR([Oll])(M0yr-^) = 1.4 x 10-^^L[oii](ergs s"^). (2) 

We correct for dust ext inction at intermedia te redshifts following 
the methods derived in Moste k et al.l ( 120121) . The methods in this 
work take into account the difference in star-formation rate between 
starburst and passive galaxies and adjusts the reddening correction 
accordingly. We then compare our SFRs with values derived from 
SED fitting (see Section l2.2.2t . We find generally good agreement. 
A straight-line fit to the relation between the two estimates of SFR 
gives a slope of 0.93, where SFRs derived from the [O II] luminos- 
ity are slightly higher than those derived from SED fitting. 

The [O II] luminosity can also be affected by the metallicity 
and gas ionisation of the host galaxy, in addition to reddening. [O II] 
is affected by the excitation state of the gas, which in turn is related 



to the ionisation (e.g. iKewlev et alj|2004 ). We do not correct for 
galaxy metallicity or the ionisation parameter, as our spectra do not 
always include [O III] or the Ha and H/? lines that are conmionly 
used to measure metallicity. In order to check that this did not affect 
the conclusions of this section, we also combine our spectra into 
bins dependent on star-formation rate derived from the SED fitting. 
Although the outflow velocities derived do vary slightly, the overall 
results of this section are not changed. We therefore continue with 
the SFRs we derive from [O II] luminosity. 

We divide our sample into four bins dependent on the derived 
star-formation rates, with the lowest bin having star-formation rates 
of < 15 M0yr~^ and the highest bin having star- formation rates of 
> 45 Moyr^^. We stack the spectra according to the method in 
Section 3.1 and fit a static and outflowing component to the Mgll 
profile using the method in Section 3.2.1. We find that the highest 
velocity outflows are hosted by the galaxies with the highest star- 
formation rates. This can be seen in Figure|8] This is consistent with 
the results already obtained, that the outflow velocity increases with 
host galaxy mass and with [O II] equivalent width (which is approx- 
imately proportional to specific star-formation rate) and suggests 
that outflow velocity depends on both parameters. It is also consis- 
tent with results from Weiner et al. (2009) who find tentative evi- 
dence for an increase in outflow velocity with star formation rate. 



10 E. J. Bradshaw et al. 



Single Gaussian 



Double Gaussian 



0.6 

1 

0.8 
O.B 

1 

0.8 
0.6 

1 

0.8 
0.6 





SFR<15M yr 



15<SFR<3CM yr"' 



30<SFR<45 M yr' 



SFR>45 M yr"' 





-3000 3000 

Velocity Offsetyi^ms"' 



-3000 3000 

Velocity Ottset/kms"^ 







Star Formation Rate vs Outflow Velocity and 


EW 


1000 






■ 


800 








. 


E 

1 600 

o 

o 400 










■ 


■5 


1 




- 








> 200 


- 


- 




■ 





■ 


0.8 






- 


-< 

1 0.6 



• - 


- 


• 
















■5 0-4 










0.2 







10 20 30 40 50 60 70 

SFR/Mgyr-' 



Figure 8. The left-hand panels show a close-up of the Mg II absoiption feature for each spectrum which have been co-added by star-formation rates, derived 
from [0 II] luminosity, along with the best-fitting single Gaussian (in red), which represents the restframe Mg II. The central panels show the breakdown of 
the absorption Une fits used to calculate the velocity of the outflowing Mgll. The red Gaussian is the non-outflowing component, the green Gaussian is the 
outflowing component, and the blue line represents the best-fit to the Mg II line (the sum of the stationary and outflowing Gaussian fits). The right panel (top) 
shows the equivalent width of the outflow as a function of star-formation rate and the right hand panel (bottom) shows outflow velocities as a function of 
star-fomiation rate. 



4.5 Specific Star-formation Rate 

We use the star-formation rates derived in Section 14.41 and the 
galaxy masses calculated in Section 2.2.2 to calculate the specific 
star-formation rate (SSFR) for the galaxy sample. We show the re- 
sults in Figure |9] A trend between excess Mg IT and SSFR is very 
apparent from the single Gaussian fits. From fitting two Gaussians, 
it is seen that the galaxies with the highest SSFRs have the highest 
velocity outflows (but again there is no trend with the equivalent 
width). 

As expected, our results suggest a strong link between intense 
bursts of star formation and the occurrence of high velocity winds. 



4.6 Stellar Population 

In this section, we use various spectral indices to determine galac- 
tic properties, and therefore divide our galaxies into sub- samples 
accordingly. First, we define the spectral indices that we will use to 
classify our spectra. 

The first of these spectral indices is the D(4000) index, which 
is a measure of the strength of the 4000A break: 



Z)(4000) 



F(4000 - 4100) 



(3) 



7^(3850 - 3950) ' 

The 4OOOA break arises from the lack of blue light from older stel- 
lar populations, the large accumulation of absorption lines from 



ionised metals, and the end of the Hydro gen Balmer series (e.g. 
iBalogh et alJI 19991 : iKauffmann et alj|2003l) . The D(4000) spectral 
index will be small for galaxies with young stellar populations and 
larger for older, metal-rich galaxies . The D(4000) index was orig- 
inally defined bv lBruzual aJ l ll983h as the ratio of the average flux 
density between the wavelengths 4050-4250 A and 37 50-3950 A . 
The definition of D(4000) from lBalogh et alj i ll 9991) that we use 
in this work is defined with narrower bands in comparison and is 
ther efore less sens i tive to reddening effects. 

ICimatti et alj 1 120081) use the mean spectral flux between 2700- 
3IOOA and 3IOO-35OOA (in the rest frame) to define the spectral 
index C(29-33): 



C(29-33) = -2.5 log 



/ f(2900) \ 
Vi^(3300)y 



(4) 



High values for the C(29 — 33) index indicate fewer very young 
stars in the galaxy, with a ga laxy classified as 'p assive' if it has 
an index (7(29 - 33) > 0.5 JCimatti et ai]|2008l) . The individual 
galaxy spectra in the UDSz with the highest (7(29 — 33) indices 
show no (or very little) [O II] emission, prominent metal absorption 
and a very red continuum . 

ISpinrad et al.l ( 119971) use a break in galaxy spectra at 2900A 
to distinguish between galaxies with older stellar populations: 
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Figure 9. The left-hand panels show a close-up of the Mg II absorption feature for each spectrum which have been co-added by specific star-formation rates, 
using the star-formation rates derived froin [O II] luminosity and galaxy masses, along with the best-fitting single Gaussian (in red), which represents the 
restframe Mg II. The central panels show the breakdown of the absorption line fits used to calculate the velocity of the outflowing Mg II. The red Gaussian 
is the non-outflowing component, the green Gaussian is the outflowing component, and the blue line represents the best-fit to the Mg II line (the sum of the 
stationary and outflowing Gaussian fits). The right panel (top) shows the equivalent width of the outflow as a function of stellar mass and the right hand panel 
(bottom) shows outflow velocities as a function of specific star-formation rate. 



The 5(2900) index is used to estimate the age of the stellar popu- 
lation where higher values of this index indicate older stellar pop- 
ulations, measuring redder colours and stronger metal absorption. 
The 5(2900) index is defined over a narrow spectral range and 
therefore it is largely independent of reddening; it is determined by 
the opacity of the m etals on the blue side of the spectral feature 
JSpinrad et al.ll 19970 . A comparison between galaxies with a high 
and low 5(2900) index is shown in FigurefTOl 

The galaxies were divided into four sub-samples for each of 
the spectral indices (in the case of the VIMOS data) and co-added 
according to the methods described in Section 3.2.1. The sub- 
samples were chosen so that they span the available range of values 
for each spectral index, but with at least 50 spectra in each stack 
in order to maximise the S/N. Data from F0RS2 were also used 
wherever possible (see appendix for details), extending the analy- 
sis towards older, redder galaxies beyond the reach of VIMOS at 
z>l. 

Figure \TT\ shows the velocities and equivalent widths for all 
three of the spectral indices defined, where the VIMOS data are 
plotted in red and the FORS2 data plotted in green. From this fig- 
ure, we see evidence for higher velocity outflows among younger 
stellar populations in both the 5(4000) and (7(29 — 33) indices. 
The trend is also apparent in the 5(2900) index, but only for the 
VIMOS sample. Unfortunately the limited number of F0RS2 spec- 
tra did not allow spectral stacks of sufficient S/N to probe the 
5(2900) index in more than two bins. 

None of the age-dependent indices show any trend with the 
equivalent width of the outflowing component in Mgll. As men- 
tioned previously, the equivalent width of Mgll is difficult to inter- 
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Figure 10. A comparison between spectra with a low B(2900) index (red) 
and the a high 52900 index (black). The difference in gradient of the con- 
tinuum around 52900 is apparent, as well as the shape and strength of the 
Mg II and Fe absorption lines. 



pret in stacked data (as the doublet is typically close to saturated) 
but this may indicate that the fraction of galaxies hosting an outflow 
is not strongly dependent on the age of the stellar populations. 
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Figure 11. This figure shows the velocity (top panels) and equivalent width (bottom panels) of the galactic outflows as a function of three spectral indices, as 
determined from two-component fits to the Mg II absorption profile in composite spectra. The green points are based on the F0RS2 sample, with the larger 
VIMOS sample shown in red. Data are plotted at the mean value of the spectral indices within each stack. 



4.7 Rest-frame U-B Colour 

We also use restframe [/ — B colours to divide the galaxies into sub- 
samples. We use a colour-absolute magnitude diagram to make a 
visual distinction between the galaxies which lie on the 'blue cloud' 
and those which lie on the 'red sequence'. We define the blue cloud 
as galaxies with U ~B < 0.5, which are generally found to be star- 
forming galaxies. Red sequence galaxies we define to have U — 
B > 0.7 and these are generally passive galaxies (although they 
can also be star-forming galaxies which are reddened by dust). The 
small population with a colour 0.7>t/ — -B>0.5is classified as 
'green valley' galaxies. 

The results of the two-component fitting analysis are sum- 
marised in Figure [12] It can be seen that the galaxy spectra show 
much the same correlation as has been seen in Section l4~6l and Sec- 
tion |4!T] The galaxies in the blue cloud show the highest velocity 
outflows, whereas those on the red sequence show the lowest ve- 
locities. 

Taken as a whole, the apparent trends of outflow velocity with 
[O II] equivalent width, spectral indices and U ~ B colour present 
very consistent trends; younger galaxies with enhanced star forma- 
tion show the highest velocity outflows. Such trends strongly sug- 
gest that outflows are driven by stellar processes (e.g. stellar w inds 
and supernovae), as discussed in lPiamond-Stanic et al.l ( l2012n . 



5 DISCUSSION 

From the outflow velocities we derive in this work, we can make a 
rough calculation of the amount of outflowing gas from the galax- 
ies in our sample. Spectra of higher resolut ion reveal that the M g II 
absorption line is typically optically thick JWeiner et al.ll2009f) . To 
obtain an estimate of the column density of the outflowing gas we 



therefore use the method outlined in ISpitzeri (Il968h . using the ra- 
tio of the Mg II doublet to estimate the optical depth. We cannot 
resolve the Mg II doublet in our data a nd hence use t he Mg II col- 
umn densities derived in this manner bv lWeiner et al.l(l2009h . as the 
average properties of blue galaxies in both samples are very similar. 

We transformed the estimates of the Mg II column density into 
a hydrogen column density using a solar abundance ratio of log 
(Mg/H) = —4.4. The Mg d epletion onto d ust grains was then ac- 
counted for (X = -1.2 dex. ISavage & Sem bach 1996) to obtain a 
hydrogen column density of N{H) = 1.3 x 10^" cm~^. 

Assuming the wind is a sphericall y symmetric thin shell of 
radius R, following IWeiner et al.l ( 120091) , it can be shown that the 
mass outflow rate will be given by the following expression: 



M ~ 44 M, 



©yr 



Ni, 



R 



(6) 



1020 cm-2 5 kpc 600 kms"^ ' 

The typical half-light radii of our galaxies are approximately 
4 kpc. A diameter of 8 kpc also corresponds to the 1 arcsec slit 
width used for our spectroscopy. Assuming the outflows are of the 
order of the galaxy radius, for an outflow of velocity of 600 km s^^ 
we obtain an outflow rate of approximately 45 M^yr^^ . 

Our estimates for the outflow rate are subject to a number of 
uncertainties, particularly since we cannot resolve the Mg II dou- 
blet. It is notable, however, that assuming a high covering fraction, 
the implied mass outflow rates in the star-forming galaxies are sim- 
ilar to the star formation rates in these galaxies (10-50 Mq yr~^), 
suggesting that outflows may play a major role in processing the 
gas content of galaxies at these redshifts. This is consistent with 
the theory that the rate of gas inflowing into a galaxy is roughly 
equal to the gas processing rate, and of this gas, half is ejected in 
the form of outflows (as discussed in this work) and half is turned 
into stars. Such behaviour is consistent with models required to re- 
produce the high-redshift mass-metallicity relation, with star for- 
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Figure 12. This figure shows the outflow properties as a function of rest-frame U — B colour The left-hand panel shows the outflow velocities determined 
from two-component fits to the composite spectra and the right-hand panel shows the equivalent widths associated with the outflows. Measurements from the 
VIMOS sample are shown in red, with the F0RS2 measurements plotted in green. Data are plotted at the mean U — B colour within each stack. 



mation an d outflows fed by the rapid accretion of external gas (e.g. 
lErbll2008l) . 

We began this work with a discussion of AGN and starburst- 
driven winds and whether either mechanism could be responsible 
for terminating star formation at high redshift. We established that 
high velocity outflows at 2 > 1 are extremely common, and driven 
by starburst activity (due to the lack of AGN signatures in our com- 
posite spectra). The outflows in many of our stacks reach velocities 
approaching 1000 km s~^ (such as those from galaxies in the high 
SSFR and high-mass samples) are likely to exceed the escape ve- 
locities of their hosts. Our findings are consistent with the work of 
Diamond-Stanic et al. (2012), who suggest that galaxies with an 
extremely high star-formation rate surface density may have sub- 
stantial momentum input from supernova and stellar winds to pro- 
duce winds of even higher velocities (>1000kms^^), without the 
need to invoke AGN-driven winds. 

In this work, we see evidence that the age of the stellar popula- 
tion correlates strongly with outflow velocity. This is shown in Sec- 
tion l4.6l where we investigate the wind properties of galaxies with 
varying spectral parameters. We find that the galaxies with younger 
stellar populations, as measured by the Z)(4000) and C(29 — 33) 
indices, as well as those galaxies with blue U — B rest-frame 
colours and high SSFR, show the highest outflow velocities. Such 
trends may be explained by a time sequence, in which all galaxies 
go through a phase of hosting high- velocity outflows following a 
major burst of star formation, which subsequently slow down once 
the galaxies evolve to have older stellar populations. 

In future work we will use morphological parameters to in- 
vestigate the role of galaxy size and orientation on the observed 
outflows (e.g. by comparing edge-on and face-on systems). Orien- 
tation effects may easily dilute the true nature of the wind velocities 
in a stacking analysis. We can also use size information to investi- 
gate influence of the star-formation surface density (see Diamond- 
Stanic et al. 2012), to investigate if more compact star-formation at 
high redshift leads naturally to more intense outflows. 



6 CONCLUSIONS 

In this work we use the spectra of 413 ii'-band selected galaxies 
to conduct a study of galactic-scale winds using a spectral stacking 
analysis at redshifts 0.71 < 2 < 1.63. We use a combination of 
data from the VIMOS and F0RS2 spectrographs, with the F0RS2 
spectra allowing us to probe older, redder galaxies at z > 1, over- 
coming a limitation of previous work. We stacked the spectra ac- 
cording to [O II] equivalent width, star formation rate, galaxy mass, 
rest-frame colours and spectral indices. We extract the velocities of 
the galactic winds by fitting a two-component model to the Mg II 
absorption profiles, representing intrinsic and outflowing absorp- 
tion. 

Evidence for outflowing material is observed in virtually all 
stacks, suggesting that strong winds may be ubiquitous for star- 
forming gala xies at these redshifts, in agreement with previous 
findings (e.g.lWeiner et al.i2009). We find that the highest velocity 
outflows are hosted by blue galaxies with young stellar populations 
and strong [O II] equivalent widths, with the highest velocity winds 
reaching velocities up to ~ 1000 km s~^. This behaviour is con- 
sistent with outflows driven by stellar processes (e.g. stellar winds 
and supemovae). The lack of AGN features in the composite spec- 
tra also suggest that AGN activity is not the primary energy source. 



Our study was motivated by the findings of iTremonti et al.l 

( 120071) ■ who observed very high velocity winds ( > 1000 km s~^) 
from highly-luminous post-starburst galaxies at 2 ~ 0.6. Extrapo- 
lating from the trends observed in our study, we speculate that such 
extreme winds may be explained as the high-luminosity tail-end 
of the distribution for normal galaxies, without the need to invoke 
feedback from an AGN. 

Assuming a high covering fraction, the implied mass outflow 
rates in the star-forming galaxies are similar to the star formation 
rates in these galaxies (10 — 50 i\fQyr^^). Such behaviour, if fed 
by the continuous accretion of cold gas, is consistent with models 
requ ired to pro duce the mass-metallicity relation at high redshift 
(e.g. lErbll2008l) . 
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APPENDIX A: FORS2 

Al Methodology 

The F0RS2 spectra have a resolution of A/ A A ~ 660 and therefore 
the Mgll doublet at 2795.5, 2802.7 can be marginally separated, 
which is not the case in the VIMOS spectra. See Section 2 in the 
main text for further information on observational details of the 
FORS2 data. 

We make a model Mg II profile by convolving two lines with a 
ratio of 1 . 1 : 1 at 2795 .5 A and 2802.7 A (to represent the absorption 
doublet) each with a Gaussian width of ct = 60kms'^. The ratio 
represents the relative depths of the two lines in the Mg II doublet 
if the absorption is optically thick, whilst the Gaussian represents 
an esti mate of the average st ellar velocity distribution, based on the 
work oflWeine r et al.l(l2009 j. We then account for the resolution of 
the FORS2 instrument by further convolving the model Mg II pro- 
file with another Gaussian with a FWHM of 5.0 A to represent the 
measured instrumental resolution, as determined from the strong 
[O II] emission in the composite spectra. 

The methodology for fitting the model Mgll profile is then 
similar to to the procedures detailed in Section 3 of the main text 
with one major difference. During the fitting procedure, approxi- 
mately 15 pixels from the centre of the Mg II profile in the compos- 
ite spectrum are excluded from the fit. This is a consequence of be- 
ing able to resolve the doublet due to the resolution of F0RS2. The 
fitting algorithm produces a more stable fit if the pixels between 
the two minima in the Mg II profile are removed. As before, how- 
ever, the non-outflowing component is fixed at zero offset velocity, 
with the depth and velocity of the outflowing component allowed 
to vary. Errors are calculated using the same methodology as the 
errors used for the VIMOS spectra; they are the result of repeating 
each fit on 100 bootstrapped spectra. 

We create a composite of all FORS2 spectra (of which there 
are 72) in order to measure the overall velocity shift of the Mg II 
profile in all galaxies and to see whether they agree with the results 
from the VIMOS data. The full wavelength range of the spectrum 
can be seen in Figure lATI The Mg II absorption from the composite 
spectrum is shown in Figure lA2l and overplotted are the best fitting 
Mg II profiles fitted using a two-component fit. The profile in green 
is the outflowing Mg II, the red profile is the intrinsic Mg II with 
its centre fixed at Okms~^. The blue profile is the sum of the two 
components and represents the overall fit to the absorption feature. 
Performing the two-component fit we measure an outflowing com- 
ponent with an offset by ~320kms~^, which is a similar to the 
velocity obtained from the VIMOS data (~340 km s~^). 
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Figure A2. This shows the Mg 11 absoiption profile from a co-added spec- 
trum (black) containing all 72 FORS2 galaxies. We use a two-component fit 
to derive the outflowing component (green) and an intrinsic component in 
red (centred at km s^ ^ ). The blue line is the overall fit, which is the sum of 
the two components. The outflowing component is centred at 320kms~^. 



We create the FORS2 composite spectra binned by the same 
parameters used for the VIMOS spectra; [O II] equivalent width, 
U — B colour, and spectral indices. The composites are not divided 
into the same sample bins as the VIMOS data, due to the differences 
in galaxy properties between the VIMOS and FORS2 samples and 
limitations in the number of galaxies per bin. For U — B colour, we 
can create a red composite at U—B > 0.7 with a significant number 
of spectra. We can create additional sub-samples for the C(29-33) 
parameter which extend to those galaxies which are extremely pas- 
sive. The extracted velocities and equivalent widths from careful 
modelling of the Mg II profiles can be seen in the main body of the 
text in Figures 4, 7 and 8 and these have been plotted against the 
results from the VIMOS data. In Figure |A3] we show an example of 
the Mg II profile fitting with FORS2 data, using U — B colour. 

This paper has been typeset from a T^U KT^ file prepared by the 
author. 



A2 Outflows in Galaxies Stacked by Property 

We divide the F0RS2 data into sub-samples and stack the data fol- 
lowing the methodology in Section [3Tl in order to investigate and 
verify the trends found using the VIMOS data. The advantage of 
using the F0RS2 data is that we can extend our analysis to those 
galaxies that are extremely red which are also at relatively high red- 
shifts. In order to measure the velocity of the outflow in each com- 
posite, we use a method similar to that in Section [3.2.1| in which 
two Mg II profiles are fitted to the Mg II absorption line. 
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Figure Al. The co-added stack of the 72 F0RS2 spectra used in this work, all selected to ensure that [O II] and Mg II were within the observed spectral range. 
Other emission and absorption features are also shown. Flux is shown in ergs s~^ cm~'^. 
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Figure A3. This figure shows the three stacks of FORS2 spectra, categorised by their U — B colour The left-hand panels show the best fitting single Gaussian, 
fit to the red side of the Mg II feature. The central column shows the breakdown of the absorption line into outflowing and non-outflowing components (green 
and red respectively), in addition to the overall best fit to the Mg II line (in blue; the sum of the stationary and outflowing components). The upper-light panel 
shows the velocity of the outflowing component for each stack, and the lower right shows the equivalent width of this component (both in green), as compared 
to the equivalent VIMOS data, which is shown in red. From top to bottom, the stacks represent the 'blue cloud', the 'green valley', and the 'red sequence' 
galaxies. 



